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Anomalous enhancements of low-energy fusion rates in plasmas: the role of ion momentum 
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Non-resonant fusion cross-sections significantly higher than corresponding theoretical predictions are ob- 
served in low-energy experiments with deuterated matrix target. Models based on thermal effects, electron 
screening, or quantum-effect dispersion relations have been proposed to explain these anomalous results: none 
of them appears to satisfactory reproduce the experiments. Velocity distributions are fundamental for the reac- 
tion rates and deviations from the Maxwellian limit could play a central role in explaining the enhancement. We 
examine two effects: an increase of the tail of the target Deuteron momentum distribution due to the Galitskii- 
Yakimets quantum uncertainty effect, which broadens the energy-momentum relation; and spatial fluctuations 
of the Debye-HUckel radius leading to an effective increase of electron screening. Either effect leads to larger 
reaction rates especially large at energies below a few keV, reducing the discrepancy between observations and 
theoretical expectations. 
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1. INTRODUCTION 

In the last ten years a number of different experiments with 
target Deuterons absorbed in a metallic matrix have found 
strong enhancements of fusion reaction rates below a few 
keV.The fusion of the implanted Deuterons with the incom- 
ing ions (D+ or Li+) has been observed; the c/(c/j/)p reaction 
has been investigated in refs. |1, 2l,[||illllllll@,[ll[nil 
and the ^'''Li(d,a)'^'^He reactions have been studied in refs. 
1(121 [Tsi [Til I15I1 finding similar strong enhancements. Ex- 
periments with gas targets show much weaker enhancements, 
which can be explained by standard electron screening with a 
potential Ue of the same order of the adiabatic limit Uad — 
28 eV. In other words, the penetration through a screened 
Coulomb barrier at energy E is equivalent to that of bare nu- 
clei at energy E + Ue- The adiabatic limit is reached when 
the correspondent enhancement of nuclear cross section can 
be explained by the gain of the electron binding energies Uad 
between the initial distant atoms and the final fused nuclei set- 
ting C/, = t/^d iQlC . 

The same type of screening mechanism could reproduce re- 
sults for deuterated metal target experiments only using an 
unreasonable large potential Ue of hundreds of eV, ten times 
greater than the adiabatic limit Uad- 

A tentative explanation ^ [H [l^ [HI], based on a sim- 
plified model of the classical quasi-free electrons, needs an 
electron screening distance of the order of the Debye length. 
This approach reproduces both the correct size of the screen- 
ing potential Ue and its dependence on the temperature: Ue °^ 
r'/^ [8, 10]. However, this model lacks a clear physical in- 
terpretation such as the one for the Debye screening, which 
is a cooperative effect, since the mean-number of quasi-free 
particles in the Debye sphere 1 17] is much smaller than one. 

The thermal motion of the target atoms is another mech- 
anism capable of increasing the reaction rate; however, 
Maxwellian momentum distribution at the experimental tem- 
peratures gives negligible effects flsl [l9[l . The relationship 



between energy and momentum of quasi particles can be 
broaden by many -body collisions |[20[l . then even a Maxwell- 
Boltzmann energy distribution leads to a momentum distri- 
bution with an enhanced power-law tail. Fusion processes 
select high-momentum particles that are able to penetrate 
the Coulomb barrier and are, therefore, extremely sensitive 
probes of the distribution tail 1,21.. , 22, ,23, ,24, 25,]. This Quan- 
tum Uncertainty Effect has been proposed in |[l9l l24l [26l [27i 
as a possible contribution to the reaction rate enhance- 
ments. 



These two mechanisms clearly cannot take into account all 
the complex many-body physics in plasmas. In different con- 
texts other approaches exist both to fluctuations and to strong- 
coupled screening il6ll . 

Among the over 50 different targets (metals and insula- 
tors) where deuteron fusion reactions have been studied, this 
present work focuses on the paradigmatic case of Ta matrix 
to show how either of the two effects we consider reduces 
the discrepancy between theoretical predictions and experi- 
mental data. This case has been extensively studied by differ- 
ent experimental groups and there exist many published data 
Jl, IS HI [^ [3; Hh in reasonable agreement. Both effects yield 
qualitatively similar results when applied to other targets, but 
detailed quantitative comparisons need further work. 

We review the experimental procedure in section |2| and 
evaluate the consequence of the Quantum Uncertainty Effect 
in section [3] The plasma screening effects on the reaction rate 
are evaluated in section |4] using for the first time the mod- 
ified Debye-Hiickel potential proposed by Quarati and Scar- 
fone in ref. [30], which is an approach introduced to study 
deviations from the weakly coupled plasma limit, where the 
standard Debye-Hiickel screening applies. We draw our con- 
clusions in section [5] 
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2. EXPERIMENTAL PROCEDURE 

All experiments performed to investigate the D + D ^ 
T + p fusion reaction, in all different target matrices, employ 
a Deuterons ion beam of energy 1 keV <Eb< 100 keV, that is 
totally absorbed by the thin foil (~ 0.1 mm) of deuterated tar- 
get, at room temperature of about 10°C. Detectors count the 
total number N{Eh,Q) of fusion reaction protons emitted in 9 
directions, and the results can be expressed in terms of the re- 
action Yield of an infinitely thick target: y^~^ — N{Ei,,Q) /Np, 
where Np is the total number of incident projectiles. 

If we consider the target nuclei at rest, the thick target reac- 
tion Yield can be expressed as: 



Y^^{Eb,Q) ^EUD — \ dE 



(2.1) 



where «£> is the density of the target deuterons, £ is the proton 
detection efficiency, o(E) is the fusion cross section and (^) 
is the energy loss for unit length, or Stopping Power. The 
effect of the thermal motion of target deuterons modifies Y^^ 
in: 



dE\ 
^1 



(2.2) 



where Vrei = |v — Vr| is the relative velocity, vi and v = 
^J2E /niD are the target and the incident particle velocities, 
E is the incident deuteron energy inside the target, mo is the 
deuteron mass, and (ov„/) — J (i^pr<I>(p()ovre/ is the thermal 
mean with ^(pf) the distribution of particle momentum tar- 
get p,. For the Stopping Power (^) all the authors adopt the 
values reported by Andersen and Ziegler 13111 . 

In the zero-temperature approximation the reac- 
tion cross section OexpiE) can be extracted from 
through a finite-interval numeric differentiation: 



, where AEh is a small 



beam energy step between two subsequent measurements of 

YZ 



.g-^p. If thermal effects are included, this procedure yields 
{{ovrei) /v)g^p instead of Oexi)- A convenient theoretical 
expression for the bare-nucleus cross section is: 



(yth{Ecm) — 



SiE,„ 



exp 




-kJ-^ , (2.3) 



where the exponential factor is the penetration function across 
the bare repulsive Coulomb potential. Eg = l/jZ^Z^e'^ / h? is 
the Gamow energy, fj is the reduced mass Zi , Z2 are the 
atomic numbers of the interacting nuclei and S{Ecm) is the as- 
trophysical factor. If the small thermal effects are neglected, 
the center of mass energy is Ecm = Et,/2. Static electron 
screening can be taken into account by introducing the elec- 
tron screening potential Ug'- 



Oth{Ec„ 



S{Ea, 



exp 



Eg 



Ue 



(2.4) 



Theoretical calculations and experimental results are of- 
ten conveniently reported in terms of the Astrophysical Fac- 
tor 5 (E): S,h,exp{Ean) = ^cm exp 



Co 



it depends much less strongly on energy. Of course, while 
Sexp{Ecm) is extracted from Y^^, the quantity Sth{Ecm) de- 
pends on the value of Ug and on the theoretical cross sec- 
tion. Sii,{Ecm) would in principle include all many-body ef- 
fects, which modify the bare cross section in vacuum. In 
this work we consider, separately, only two effects as dom- 
inant for the reaction rate: the quantum broadening effect 
on the momentum distribution and the electron screening, 
which we evaluate beyond the Debye-Hiickel approximation 
given the plasma density inside the metal matrix. The case of 
Deuterated Tantalum target has been investigated extensively 
m 3, m 1^ I7, 8]. Measurements are reported with great de- 
tail in ref. |6], then we adopt the relative experimental en- 
vironment for our comparison: matrix target chemical com- 
position is pure Ta (Z = 73, A= 180.948 a.m.u.), target tem- 
perature T = 10°C, Ta density at room temperature p=16.65 
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g cm (riTa = 5.54- 10 cm ), stoichiometric coefficient 



: 7.9, absorbed Deuterium density no 



0.701 -1022 cm 
relative to the 

deuterated Tantalum experiment of ref. |6]. The experimen- 
tal results can be reproduced only with an electron screening 
potential Ug — 309 eV ^ Uad — 28 eV; similar large screen- 
ing potentials are needed to fit other deuterated-Tantalum ex- 
periments: Ug = 322 eV in ref. |3] and Ug = 340 eV in ref. 
yTlla]. The experimental procedure has been critically studied 
in 0, Btl, especially target surface contaminations and inho- 
mogeneities in the implanted deuterons distribution have been 
identified as possible sources of systematic errors in the ex- 
perimental data. However, performing deuterated Tantalum 
targets experiments with different surface contaminations, the 
electron screening potential needed to reproduce the data is 
in the range Ug = 210 — 460 eV much greater than 

Uad and in agreement with the previous results. In the ab- 
sence of incoming beam there are ions and conduction elec- 
trons inside the metallic matrix, while the ion beam may ion- 
ize deuteron atom target or metallic ions, so additional charges 
may appear Accurate numerical simulations performed by 
Huke et al. |32] suggest migration of electrons from the host 
metal atoms to the Hydrogen ions during the impact event. 
This effect produces a screening potential Ug = 39.7 eV for 
deuterated-Tantalum target, greater than Uad but still one or- 
der of magnitude lower than the one needed to explain the 
experimental results. 



Fig. [U shows both Y^^ and the extracted Sgxp, 



smce 



THERMAL EFFECTS AND QUANTUM 
UNCERTAINTY. 



Collision frequency among target and incoming deuterons 
is the most important quantity to consider in this treatment. 
Collision frequencies among other particles, e.g., target atoms 
and metalUc ions, are less important and are not responsible of 
modification of momentum distribution function and of fusion 
rate |33]. Momentum-energy dispersion relation is broad- 
ened by Quantum Uncertainty Effect (QE), so that a Maxwell- 
Boltzmann energy distribution 'PiE) oc exp{—E/ki,T) can 
yield corresponding momentum distributions ^(p;) with en- 
hanced tails. Since nuclear fusion cross section depends on 
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FIG. 1: Experimental results from ref. f3|. Upper panel: Y^^ data points 
from one of the 13 runs weighted to extract the cross section Ogyp ; data nor- 
malized to a beam energy Ej, (Lab frame) of 50 keV; curves from eq. 
i2.l\ are computed with eq. (23) (bare) and eq. i2A\ (screened) cross sec- 
tions. Lower panel: Astrophysical Factor data points S^xp extracted from Y^p ; 
bare Astrophysical Factor is assumed to be 5(£„ii) =43+ 0.54£„ii keV b, 
by normalization with higher energy results (£„„ > 40 keV), while screened 
5,1, curves are obtained from eq. i2.4i . Thermal effect are neglected. 



the relative momentum e^^^, = j/jvve/ such an effect increases 
the value of {<3{Zp^.^i)vrei) /v respect to the one obtained with 
a sharp energy-momentum relation. It is a good approxima- 
tion of the experimental situation to consider a beam of par- 
ticles with definite energy and momentum and a Maxwell- 
Boltzmann energy distribution of thermalizedmtarget parti- 
cles. The momentum-energy relation of the target particle can 
be represented with a Lorentzian f^^, at least for E, ^ Ep-. 

where e^, = ^niovj , y= hnofScoUVcoii , ^coU and Vcou = 
\/2Etlmp are the collisional cross section and velocity. In 
ifioi 13, IS llSl Quantum Uncertainty Effect has been ap- 
plied to the fusion reaction between beam and absorbed 
deuterons, adopting a Coulombian collisional cross section 
'^coU ~ ^"^/^p,- This QE produces a rate increase, but only 
at beam energies lower than the energy at which the experi- 
mental results start rising {Eb ^ 2 keV for ref. [6] instead of 
Eh ^ 6 — % keV, as shown in fig. |2]l. 

Zubarev ll28l I29I1 proposes that the QE should be effec- 
tive only for a small fraction of target deuterons that are in 
a quasi-free mobile plasma states in the reaction zone and that 
the rate enhancement is due to reactions of deuterons of this 
small plasma fraction with the other stationary deuterons in 
the target. The beam is only needed to maintain the plasma 
states. Zubarev reproduces the enhancement with this mecha- 
nism and a Coulombian collisional cross section. 

However, collisional cross sections depend on the plasma 
environment. In particular, the Coulombian cross section 

<^coU = l^p, is appropriate only for weakly interacting plas- 

2 

mas, i.e., plasmas with parameter F = j^fj— <C 1, where 

aj„. = (3/47tn£))'/^ is the Wigner-Seitz radius, but T 100 
for the typical experimental conditions we are interested in, 
see for instance ref. [6], where no — 0.701 • 10^^ cm^^, 
khT = 0.0244 eV, = 3.24 • 10"^ cm, resulting in T = 182 
(very close to the liquid-solid transition). Therefore, a strong 
coupled plasma scheme should be used for the absorbed 
deuterons. A model that can describe the screened colhsional 
cross section in this limit is the Ion Sphere Model ITtII : 

a,„zz = 27t a? fl2,. (3.2) 

Eq. 13.21 is a non trivial expression for QcoU, because the pa- 
rameter ai contains information on the particle-particle cor- 
relations and in ref.s Il33.i35ll (and citation therein) Ichimaru 
and collaborators have given a detailed analysis of its expres- 
sion and evaluation. For what concerns the present problem 
we can say that ai is a correlation factor of the order of unity 
or, more precisely, that its value lies between 0.4 and 0.9. 

The resulting distribution for p, is 't'(Pf) = 

^JodErbyiE, - E,„)e-^'/'"'T, where In = 

4% pjdpt dEt?>y{Et — Ep,)e^^'l'^i'T jg (he normahzation 
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integral. Therefore, the thermal mean is: 



300 



' iicosi> 



(3.3) 



where ■& is the angle between p, and the beam and v^/ — 
p, cos-d)'/ . 

Eq. ( 13.31 1 has been numerically evaluated with parameters 
appropriate to the experiment of ref. 16|] using a constant As- 
trophysical factor 5(6;,,.^, ) ~ 5o = 43 keV b, valid within w 6% 
for E[, < IQ keV. The ■& integration has been done analyti- 
cally in terms of incomplete Gamma-Euler function, while the 
remaining integrations have been performed using the Gauss 
adaptive method; results are shown in fig. |2]for Coulombian 
and Ion Sphere Model (ISM) collisional cross sections. Both 
cases yield a strong enhancement of the reaction rates at low- 
energy: below Eh ^ 2 keV for the Coulombian and below 
~ 10 keV for the ISM case. These behaviors should be 
compared to the anomalous enhancement of the experimen- 
tal data that start below energies ~ 4 — 6 keV. In the ISM 
case, the energy-threshold below which the rate is enhanced 

-2/3 

depends on Ocou that, in turns, depends on ' , (see eq. 
( I3.2l i). Our calculation uses no adjustable parameter (ai — 1 
in fig. O, however, if we follow the approach of Kim and 
Zubarev in ref. |29J, we could apply the enhancement only to 
that fraction of the absorbed deuterons that is in a quasi-free 
plasma state. The dependence of the reaction rate r on the 
fraction / = no,,, /no of absorbed target Deuterons, no, that 
are quasi-free, no,,,, is not trivial: r is not simply proportional 
to /, but it has an additional dependence on / through the 
ISM collisional cross section Ocoii eq. ( 13.21 ); it is clear, how- 
ever, that / = 1 gives the maximal effect. In principle, the 
fraction / could be calculated in some theoretical model, de- 
termined from an independent experimental measurement, or 
used as a free parameter to be fitted. Results presented have 
been obtained with f —I and, therefore, should be interpreted 
as upper limits for this effect. 



4. MODIFIED DEBYE-HUCKEL SCREENING. 

We discuss the effect of electron screening on deuteron- 
fusion reactions in metal matrix. We do not consider the in- 
terplay with the effect discussed in the previous section. The 
screening potential at distances near and below the turning 
point influence the most reaction rates. Thermal effects and 
screening phenomena in plasma environment are strictly con- 
nected. Assuming, in a ideal plasma, a Maxwell-Boltzmann 
(M-B) thermal distribution, the Debye-Hiickel ion-screening 

potential VoHir) = ^'^^'^ exp{—r/RoH) is obtained by means 
of the linear Poisson equation (Ron = ki,T / A%nZe^ is the 
Debye-Hiicke radius). In strongly coupled plasmas, Debye 
description loses its physical interpretation. In the case stud- 
ied, the main reason is that the number of particles inside the 
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FIG. 2: Astrophysical Factor experimental points from ref. f^; is the in- 
cident Deuterons energy. Continuous curves sliow tlie zero-temperature As- 
tropliysical factor S for screened (tliin curves) and bare (thick curves) nuclei; 
here we adopted for S the same expressions of fig. [T]with = Et/2. The 
dashed and dot-dashed thick curves represent the thermal mean with Quantum 

Uncertainty Effect 5 = {ov„i)qe ^ y''^exp ^Jty'^^ j in keV b, given by 
the numerical integration of eq. )3.3t : specifically, the dot-dashed curve is rel- 
ative to the Coulombian c,-,,;;, while the dashed one is plotted for the Ion 
Sphere Model Ccoii, e<5- 13. 2K with ai = 1. 



Debye sphere is too small. The D-H approach needs to be 
extended to include features that appears as systems leave the 
weak-coupled regime. Near the weak-coupled regime devia- 
tions from the D-H regime can be parameterized by an elec- 
tron cloud with a steady-state generalized spatial distribution 
of q-type (for q ^ I electrons are distributed according to a 
Boltzmann factor): we assume that such deviation can be an- 
alytically continued in the strong-coupled regime. Following 
this strategy, Quarati and Scarfone 113011 have recently derived 
a new screening potential called Modified Debye-Hiickel po- 
tential using two different approaches. The first one uses a 
generalized non-linear Poisson or BernuUi equation, the sec- 
ond is based on superstatistics ll36l[37ll . We discuss this sec- 
ond approach in some detail in the following. 
We choose the value ^ = 0, also for the Ta matrix, because we 
can reproduce an electron distribution spatially concentrated 
around the deuteron, with a strongly depleted tail with cut-off 
at Roh/^- We are encouraged in this line of treatment also by 
the result of ref. llioll . 

The authors of ref. jsoll assumed that non-linear ef- 
fects produce fluctuations on the inverse Debye-Hiicke ra- 
dius 1 /RoH, with a Gamma-function probability distribution: 

1 



.fq{rXh)) = 



where/^(r,A.,A,o) 



represents the probability density to observe a certain value A- 
spreads around a central value Xq. To obtain from fq{r,X,Xo) 
an electron depleted tail distribution with a cutoff, we limit the 
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entropic index q into the < <7 < 1 interval, assuming: 



1 

(1 -?)^(?)!*^o 



where g{q) is a generic entropic index function, that satisfies 
the condition g(l) — 1 (inref. 13^ the choice g(g') — {2~q)^^ 
is adopted). The point charge potential y^(r) can be identified 
by the functional: 

through the relation: 



D 



Ca 



1 



(4.1) 



where C, 

„2 



normalization factor, D 

I 1 



- (2 - q)^{q) i 

ZiZ^e"-!^-^ and Xq = (A.) = /o°°/,(r, A.,Xo)A.<iA. - 
The charged particles distribution p (r) and the point charge 
potential Yq{r) °^ p('") can be derived from eq. ( 14. Il l, develop- 
ing the functional y n(rA n) with the previous assumption. A 
Tsallis cut-off form Il38l[39ll is obtained for the potential: 



Z1Z2 



if r < 



(4.2) 



if r> 



We note the relation of the power-law expression of Yq{f) 
with Tsallis-like distributions. In systems with long-range in- 
teractions and correlations and/or memory or in systems that 
are not in global thermodynamic equilibrium, but rather in a 
metastable state with a long, finite lifetime, such deformed 
Tsallis-like distributions are relevant at the very least as con- 
venient parameterizations of deviations from equilibrium; su- 
perstatistics is an example of such an approach ll36[|37ll . 
The main contribution to the charged particles fusion cross 
section is given by the screening barrier penetration factor 
P(E). In the standard Debye-Hiickel potential case, the sim- 
plified expression o(E) = ^^P{E) can be obtained, that dif- 
fers from the bare nuclei cross section of eq. i2.3i only 

for the penetration factor P{E) — exp 
Udh 



Eg 
E+Udh 



where 



Rdh 



. In the case of the modified Debye-Hiickel po- 
tential Vq{r), the penetration factor Pq{E) is given by the ex- 
pression: 



Pc,{E) = exp 



4f [2..^(V,M-£)]i 



dr 



(4.3) 



where the classical turning point has to be determined 
through the equation Vq{rq) — E. 

Eq. (14.3b can be analytically solved in the q = Q 
case, obtaining rq = ZiZ2e^ / [E + g{Q)D) and Po{E) = 



exp I 



E+g(Q)D ) ■ 
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FIG. 3: Astrophysical Factor experimental points from ref. f^. Bare nuclei 
curve correspond to Stare = 43 + 0.54£'„ii keV b, while the screened curves 
are5 = /,-5'4„„(£„n) , where is defined in eq.g3), J/, = e- (^) ,g{q) = 
i — 2q and ^ = 0. Ion thermal motion is neglected: £„ji = Ei,/2. 



The bare nuclei cross section reported in eq. ( 12.31 ) can be cor- 
rected, to account the modified Debye-Hiickel (D-H) screen- 
ing, multiplying <3hare{E) by the factor 



fq = Pq{E) exp %\ — 




(4.4) 



For instance we compare the Modified D-H Astrophysical 
Factor for Dh-D reaction with the ref. {§\ experimental data, 
adopting the choice g{q) = 3 — 2q. The results is shown in fig. 
[3]for the entropic index ^ = 0. One can observe as a screening 
potential Uq — D, three times lower than in the standard D-H 
case (Uq ^ 100 eV instead of Udh ^ 300 eV), is required to 
reproduce the experimental data. An electrostatic screening 
potential of this order of magnitude has been obtained, for in- 
stance, by Saltzmann and Hass |40] through a Thomas-Fermi 
model of the electron gas in a deuterated-copper target (they 
obtained a screening potential of 163 eV, instead of the 470 
eV needed to reproduce the experimental results). 
Now we are treating g{q) and q as free parameters, but, in 
principle, a link can be establish between the inverse D-H Ra- 

dius, the temperature fluctuations and the <7-index: ^ ^ = 

Rdh 

^f^Y^ — Vl ~ 1 by the electron charge ^-distribution 
around the ion (if experimentally known or as deduced from 
a model). By this way the modified Debye-Hiickel potential 
can be obtained starting from the environment condition. 



5. CONCLUSION 

We studied two effects of non-Maxwellian velocity dis- 
tributions in plasmas as possible explanations for the fusion 
rate enhancements observed in deuterated-metal target exper- 
iments: modifications to the thermal mean (ovve/) due to the 



6 



Quantum Uncertainty Effect and the stronger screening aris- 
ing from non-thermal electron distributions. Either effect re- 
duces the discrepancy between theoretical models and exper- 
imental results. 

The broadening of the momentum distribution due to the 
QE has been studies using a collisional cross section Ocoii de- 
rived by the Ion Sphere Model, eq. (13.2b . model that captures 
the main features of collisions in the strong coupled plasmas 
that characterize the experimental environment. 

We have been able to obtain a rate enhancement up to en- 
ergy thresholds much higher than previous calculations (zt!] 
and even above the one observed in experiments, as shown in 
fig. 121 If one makes the reasonable additional hypothesis that 
only a fraction of the absorbed deuterons are in a quasi-free 
plasma-state, which can be shown to be still strongly coupled, 
one should be able to reproduce the experimental behavior 



We intend to explore more quantitatively this possibility in 
the near future. 

We have also considered the effects of plasma screening on 
the reaction cross sections, when the modified Debye-Hiickel 
potential introduced by Quarati and Scarfone in ref. [3^] is 
used. Already in ref. ||40.1 a screening potential Ue ^ 100 eV 
was obtained, greater than the adiabatic limit Uad, but still too 
low to reproduce the experimental results. In this work we 
showed that the modified D-H potential of eq. ( 14.21 ). with an 
appropriate choice of the function g{q) and of the entropic in- 
dex q, can reproduce the data, see in fig. [3] In principle the 
value of the entropic index q can be derived from the plasma 
equation of state and g{q) from the electron charge distribu- 
tion around the ion, then the modified D-H screening contribu- 
tion to the rate enhancement can be evaluated without free pa- 
rameters. We shall investigate this last point in a future work. 
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